Rydberg-induced Solitons: Three-dimensional Self-trapping of Matter Waves 



O 

(N 

u _ 

Ph. 
< 

(N 



Oh' 
I 



O 

Oh' 



(N 
O 




F. Maucher/ N. Henkel/ M. Saffmaii,^ W. Krolikowski,^ S. Skupin,!'"* and T. Pohli 

^ Max Planck Institute for the Physics of Complex Systems, 01187 Dresden, Germany 
^Department of Physics, University of Wisconsin, Madison, Wisconsin 53706, USA 
^ Laser Physics Centre, Research School of Physics and Engineering, 
Australian National University, Canberra, ACT 0200, Australia 
^Friednch Schiller University, Institute of Condensed Matter Theory and Optics, 077^3 Jena, Cermany 

We propose a scheme for the creation of stable three-dimensional bright solitons in Bose-Einstein 
condensates, i.e., the matter-wave analog of so-called spatio-temporal "light bullets". Off-resonant 
dressing to Rydberg n_D-states is shown to provide nonlocal attractive interactions, leading to self- 
trapping of mesoscopic atomic clouds by a collective excitation of a Rydberg atom pair. We present 
detailed potential calculations, and demonstrate the existence of stable solitons under realistic ex- 
perimental conditions by means of numerical simulations. 
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Self-trapped nonlinear waves and the possibility to cre- 
ate "particle-like" wave packets have fascinated scientists 
over the last decades [1-4]. In nonlinear optics, the 
creation of stable three-dimensional bright solitons — so 
called "light-bullets" [5] — has been under active pur- 
suit [4], but was realized only recently in a discrete set- 
ting of waveguide arrays [6]. One major obstacle stems 
from the fact that nonlinear confinement usually comes 
hand in hand with collapse instabilities [7]. In prin- 
ciple, this problem can be overcome via nonlocal non- 
linearities, where the nonlinear self- induced potential at 
a particular point in space depends also on the nearby 
wave amplitudes. While experiments have established 
two-dimensional optical solitons in a continuous nonlo- 
cal medium (see e.g. [8]), the realization of their three- 
dimensional counterparts remains an elusive goal. On 
the other hand, Bose-Einstein condensates (BECs) have 
emerged as clean model systems to investigate soliton 
formation [3, 9, 10], and, more recently, to study the 
effects of nonlocal nonlinearities, as arising from dipolar 
interactions [11]. Such interactions can give rise to stable 
two-dimensional bright solitons [12], while the partially 
attractive nature of dipolar interactions precludes stabil- 
ity of three-dimensional solitons due to collapse [13]. 

In this work, we propose an experimental scheme for 
the creation of "matter wave bullets", i.e., stable three- 
dimensional (3D) bright solitons in a Bose-Einstein con- 
densate. The approach is based on optical dressing 
of a groundstate atom BEG to highly excited Rydberg 
states [14-16] (see Fig. 1). Within detailed potential cal- 
culations for Rubidium condensates we identify an appro- 
priate range of Rydberg states that provides attractive 
nonlocal nonlinearities, which do not lead to condensate 
collapse. We present extensive 3D simulations of the un- 
derlying Gross-Pitaevskii equation that demonstrate the 
emergence of stable condensate bullets for realistic ex- 
perimental parameters and scenarios of their creation. 
The calculations show that mesoscopic numbers of ul- 
tracold atoms can be bound together by a single pair of 



FIG. 1. (a) Simplified level scheme of the considered Rydberg- 
dressing approach. The off-resonant coupling of condensed 
ground state (\g)') atoms to strongly interacting Rydberg 
states (|e)) gives rise to a stable bright soliton bound by a 
rim of a collective quantum excitation of a Rydberg atom 
pair, shown in (b). Depicted is a numerically obtained soli- 
ton formed by dressing a Rb condensate to 65-D3/2 Rydberg 
states with a Rabi frequency Sl/27r = 0.5MHz and laser de- 
tuning ^/2ti = 32MHz. Shown are half-maximum iso-density 
surfaces of the BEG (inner red sphere) and the Rydberg-pair 
density (outer blue shell), while the gray-scale gives the inte- 
rior density. 



strongly interacting Rydberg atoms which is coherently 
shared throughout the condensate (see Fig. 1). This ex- 
otic state — arising from intricate coupling between in- 
ternal and translational atomic quantum dynamics — is 
shown to live over hundreds of milliseconds and could, 
thus, provide the first realization of 3D bright solitons in 
atomic BECs. 

In the zero temperature limit the BEG wavefunction 
■0 obeys the Gross-Pitaevskii equation (GPE) 
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where N — {tp\tp) is the number of atoms, M denotes the 
atomic mass and g — Anh^a/M describes the strength of 
local interactions due to s-wave scattering with scattering 
length a. The range of the potential W greatly exceeds 
the average interatomic distance, and, thus gives rise to 
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the nonlocal nonlinearity in Eq. (1). Let us first consider 
a simplified, spherically symmetric potential of the form 
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Typical values of the parameters Ce < and i?c will 
be discussed below in our explicit calculations of the 
dressing-induced interaction. As we will see, this generic 
potential readily captures the essential physics of soliton 
formation. Introducing appropriate length (i?c) and time 
(r = RlM/h) scales, Eq. (1) can be re-expressed as 
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This simple scaling shows that the BEC properties are 
solely determined by two dimensionless parameters a = 
-MNCq/H^RI and 7 = AiraN/Rc, where {tpltp) = 1. 
In order to find approximate soliton solutions we first 
use a Gaussian trial wavefunction of width a. Following 
standard variational analysis we find an implicit relation 
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involving the interaction ratio r/ — 7/0; and the width 
of condensate a. The integral Z — —J W{r — 
r')e-'-'/'^'e-'-"/'^'dr'dr can be evaluated analytically in 
terms of hypergeometric functions. In Fig. 2(a) we show 
the resulting existence curves for different values of the 
interaction ratio rj. As expected, the condensate size is on 
the order of and increases with increasing local repul- 
sion. According to the Vakhitov-Kolokolov criterion [18], 
the stability of the solitons is expected to change at the 
minima of the existence curves. This is confirmed by 
our numerical solutions of Eq. (3) , which show that soli- 
tons on the left branches are stable, whereas the right 
branches are unstable. Consequently, the minima mark 
the critical a above which stable self-trapped solitons can 
exist. The corresponding stability diagram is shown in 
Fig. 2(b). For vanishing local repulsion (7 — 0) this crit- 
ical value is found to be a = 5.7, and it increases almost 
linearly with 7 in the considered parameter range. 

While this analysis yields a simple picture for the sta- 
bility conditions of self-trapped BECs, their dynamical 
formation deserves some additional discussion. In par- 
ticular, we note that attractive kernels generally lead to 
long wavelength modulation instabilities, i.e an exponen- 
tial growth of periodic perturbations on a constant den- 
sity background. The wavenumber fcmi of the periodic 
perturbation with maximum growth-rate gives the typi- 
cal length scale Ami — 27rfc~j^ of the unstable modulation. 
As shown in Fig. 2(a), the soliton width a of the stable 
branch is significantly smaller than fcmi- Hence, there is 
a wide range of initial condensate states with diameters 
smaller than Ami and with sufficiently large values of a 
to eventually form a stable soliton. 




FIG. 2. (a) Existence curves for stable (thick solid lines) and 
unstable (dotted lines) solitons for several rj = y/a and typical 
length scales of modulational instability (thin blue lines) (b) 
Corresponding parameter regions for self-trapping. The dots 
show numerical results using the potential Eq. (5) for a Rb 
condensate with atom numbers A'^ — 1000 , 3000 , 4850. 



Having established the general conditions for self- 
trapping based on the simplified potential Eq. (2), we 
now discuss its physical realization via off-resonant dress- 
ing of ground state atoms with high-lying Rydberg states. 
Let us consider first the most simple case of a sin- 
gle non-degenerate Rydberg state coupled to the atomic 
groundstate with a Rabi frequency ^l/2n and laser de- 
tuning A/27r [see Fig. 1(a)]. This admixes a small frac- 
tion V — (i7/2A)^ of Rydberg character into the atomic 
ground states. For two distant atoms, this leads to an ef- 
fective interaction v^C^jr^ ^ arising from the strong van 
der Waals interaction C^/r^ between Rydberg atoms. At 
smaller distances r the doubly excited state is blocked 
by the van der Waals shift [17], such that the effective 
interaction saturates. Altogether, this yields an effec- 
tive potential of the type Eq. (2), where Cq — v^Cf, and 
i?e - (-C6/2ftA)i/6 ]14]. 

This scenario may be realized with nS 1/2 Rydberg 
states of Rubidium atoms. However, the van der Waals 
interactions of alkaline atoms are entirely repulsive for 
highly excited S'-states, such that the above scheme 
only produces repulsive nonlinearities. We, therefore, re- 
sort to higher angular momentum {nPj or nDj) states, 
accessible either by single- or two-photon transitions. 
The underlying Rydberg-Rydberg atom interactions are 
calculated by diagonalizing the two-atom Hamiltonian, 
which includes the dipole-dipole coupling between Ryd- 
berg states. Since there are {2J+ 1) Zeeman degenerate 
Rydberg states {mj = — J, — J + 1, J) one obtains a 
set of (2 J + 1)2 potential curves vjf^{r) [see Fig. 3(a)] 
with respective molecular eigenstates |^^(r)) ]19]. For 
distances r larger than the van der Waals radius i?vdw 
]19] the potentials behave like v}y^J{r) - V^^- The 
eigenstates depend on the interatomic distance r and on 
the angle between the molecular axis r and the quan- 
tization axis defined by the coupling laser. Note that all 
of the potential curves V^^^ (r) must be of negative sign. 
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3. (a) Two atom level-scheme for dressing to 65-D3/2 

Ryd, 



Rydberg states and Rydberg interactions V^'^ (r). (b) Re- 
sulting single potential curve {rhj — 3/2) for 1} — according 
to Eq. (5). (c) Attractive and repulsive branches of interaction 
strengths Cg^' vs. n. For n > 59, all curves are attractive. 



Already a single repulsive curve would lead to resonant, 
simultaneous excitation of two Rydberg atoms at a dis- 
tance for which v}^^{r) = 2hA, and cause rapid loss on 
a timescale ^ 100/is for typical Rydberg states (see be- 
low). This requirement significantly restricts the range 
of appropriate principal quantum numbers n and angu- 
lar momentum states. We have performed potential cal- 
culations for different Rydberg states of Rubidium, and 
found n£'3/2"iiiteractions to be most favorable, as they 
yield entirely attractive curves for n > 59 [cf. Fig. 3(c)]. 

The corresponding two-atom level scheme for dress- 
ing to such Rydberg states is shown in Fig. 3(a). Upon 
proper choice of laser polarization, the atomic ground 
state {\g) = |5S'i/2)) can be selectively coupled to a sin- 
gle mj-state out of the degenerate nD3/2-manifold [cf. 
Fig. 4(a,b)], denoted by |e) = \nD3/2{mj)). Generally, 
the singly excited two atom states (|ge), \eg)) are, how- 
ever, coupled to all molecular eigenstates |/Li/3(r)), due to 
TO-mixing by the dipole-dipole interaction D{r). The cor- 
responding Rabi frequencies 17^ = {ee\fj,p{r))fl are given 
by the overlap with the laser-coupled Rydberg states, and 
thus depend on the interatomic separation vector. 

Assuming il <C A the corresponding two-atom Hamil- 
tonian can be diagonalized within fourth order perturba- 
tion theory. Omitting terms that do not depend on r, one 
obtains a single anisotropic effective interaction potential 
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between two dressed ground state atoms. This two-atom 
calculation can be straightforwardly extended to the A^- 
atom case. As long as the total number of excited Ryd- 
berg atoms Nryd = (^)^-^ < 1 [20], this simply yields 
J2i<j^i\^i ~ for the total potential energy of A^ 
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FIG. 4. Effective ground state potentials (c,d) resulting from 
two- photon state dressing (a,b) to rhj = 3/2 (a,c) and 

rhj = 1/2 (b,d) Rydberg states, with a two-photon Rabi fre- 
quency Q./2ti — {QiQ2/2S)/2n = 0.5MHz and laser detuning 
A/2tt = 32MHz. The aspect ratio of the resulting soliton is 
shown in (e) as a function of for rhj = 3/2 (circles) and 
rhj = 1/2 (squares). The insets in (e) show respective 3D soli- 
ton profiles [cf. Fig. l(b)[ for ft/2TT = 0.3MHz and 0.5MHz. 



atoms at positions [14]. The mean field dynamics of 
the BEG is, thus, governed by the GPE (1) and Eq. (5). 

To further illuminate the physical mechanism behind 
the observed self-trapping, let us consider the resulting 
many-body states of the Rydberg-dressed condensate. 
From the described perturbation theory we find that the 
second order contribution to the many-body groundstate 
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is a coherent superposition of Rydberg atom pair states, 
whose correlation function drops to zero within a dis- 
tance for which V^^^^{r) — —2hA. This radius coincides 
with the typical size of the self-trapped BEG, such that 
two atoms are excited on opposite sides of the soliton. 
According to Eq. (6), these pairs are coherently shared 
among all atoms, resulting in a collective shell of Ryd- 
berg excitations ]c.f. Fig. 1(b)]. It is the strong attraction 
between these Rydberg atom pairs that self-confines the 
BEG. 

In the following, we present explicit calculations for a 
^^Rb condensate of A^ = 1000 atoms dressed to 650^/2 
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FIG. 5. Self-trapped condensate dynamics after instantaneous 
turn-on of the dressing lasers and sudden release from a har- 
monic trap (wtr = 170Hz). The white lines give the width of 
a freely expanding BEG. 651)3/2 Rydberg states are coupled 
with n/27r = 0.5MHZ and A/27r = 32MHz. 

Rydberg states with A/27r = 32MHz, based on a 3D in- 
tegration of Eqs. (1) and (5). To ensure applicability of 
our potential calculations, we consider Rabi frequencies 
between IMHz and 3MHz for which 0.006 < A^^yd < 0.06. 
Figs.4(c,d) show that due to the angular dependence of 
the underlying Rydberg wave functions the resulting in- 
teraction potentials are anisotropic [in contrast to our 
isotropic model potential Eq. (2)], with a cylindrical 
symmetry around the quantization axis (z-axis). Con- 
sequently, the resulting self-trapped ground states shown 
in Fig. 4(e) are also slightly asymmetric. However, the ki- 
netic pressure as well as the local contact interaction both 
tend to round out the soliton solutions, such that our 
previous analysis using the simplified potential Eq. (2) 
is expected to remain valid. In order to make quanti- 
tative comparisons, we define an angular dependent ra- 
dius Rci'd) at which the potential Eq. (5) assumes half 
its minimum value W{0) — —hil'^/8A^. For the cho- 
sen parameters this radius is considerable larger than 
Rvdw ~ 3/im, such that the actual potential Eq.(5) 
closely resembles our model potential Eq. (2). Using the 
average softcore radius = tt^^ Rc{'&) sin ddi} as a 
natural length scale and f — R^M/h to scale times, we 
obtain a dimensionless GPE of the type of Eq. (3), with 
a = Nn^R'^M/8hA^ and 7 = AnaN/Rc. The numer- 
ically obtained critical interaction strengths for soliton 
existence, shown in Fig. 2(b) are in excellent agreement 
with the variational results based on the simple model 
potential Eq. (2). 

Experimentally, coherent Rydberg excitation of cold 
atoms has been demonstrated under various conditions 
[21]. Ideally, the preparation would start from a cylin- 
drically symmetric trap to initialize a BEC already close 
to the solitary ground state. However, in order to demon- 
strate robustness of the soliton. Fig. 5 shows the dy- 
namics of a significantly different initial condensate af- 
ter sudden release from a spherical harmonic trap and 
simultaneous turn-on of the dressing lasers. Compared 
to the undressed BEC expansion, self-trapping becomes 
evident within 10 ms, manifested in stable oscillations of 
the condensate size. On the other hand, the lifetime of 



the dressed BEC is as long as several 100 ms, being pre- 
dominantly limited by the decay of the weakly admixed 
6-P1/2 state [14] ]cf. Figs.4(a,b)]. A verification of the pre- 
dicted self-trapped "matter wave bullets", thus, appears 
to be well in reach of current experimental capabilities. 

In conclusion, we have shown that off-resonant dressing 
of BECs to attractively interacting Rydberg states pro- 
vides a promising route for the first realization of stable 
self-trapped three-dimensional bright solitons. While we 
chose Wa{nD^/2) atoms as one relevant example, the pro- 
posed scheme generally applies to Rydberg states of any 
atomic species with sign-definite attractive van der Waals 
interactions. Explicit potential calculations and numeri- 
cal simulations of the resulting GPE have been performed 
to demonstrate experimental feasibility and to work out 
appropriate parameters. In addition, we have shown that 
a simple isotropic model potential captures the essential 
physics of the observed soliton formation, which may be 
useful for future theoretical studies on, e.g., 3D soliton 
interaction or higher-order states [22]. The anisotropy 
of the interaction together with its tunability may also 
open up new routes to transfer angular momentum to the 
condensate ]23]. 
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